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We present molecular dynamics simulations of driven vortices in layered superconductors in the presence of 
an external homogeneous force and point disorder. We use a model introduced by J.R.Clem for describing 3D 
vortex lines as stacks of 2D pancake vortices where only magnetic interactions are considered and the Josephson 
interlayer coupling is neglected. We numerically evaluate the long-range magnetic interaction between pancake 
vortices exactly. We analyze the vortex correlation along the field direction on (c-axis). We find that above 
the critical current, in the "plastic flow" regime, pancakes are completely uncorrelated in the c-direction. When 
increasing the current, there is an onset of correlation along the c-axis at the transition from plastic flow to a 
moving smectic phase. This transition coincides with the peak in the differential resistance. 



1. INTRODUCTION 

The prediction 10 of a dynamical phase transi- 
tion upon increasing drive, from a fluidlike plastic 
flow regime [0 to a coherently moving solid § in 
a moving vortex lattice has motivated many re- 
cent theoretical [0, experimental |4|-|(| , and sim- 
ulation work In a previous work we have 
studied the dynamical regimes in the velocity- 
force curve (voltage-current) in 2D thin films |L0| 
and found two dynamical phase transitions above 
the critical force. The first transition, from a plas- 
tic flow regime to a smectic flow regime, is charac- 
terized by the simultaneous occurrence of a peak 
in differential resistance, isotropic low frequency 
voltage noise and maximum transverse diffusion. 
The second transition, from a smectic flow regime 
to a frozen transverse solid, is a freezing transition 
in the transverse direction where transverse dif- 
fusion vanishes abruptly and the Hall noise drops 
many orders of magnitude. In other 2D simula- 
tions the peak in differential resistance was found 
to coincide with the onset of oricntational order 
D and a maximum number of defects J^] . Exper- 



imentally, the position of this peak was taken by 
Hellerqvist et al. j^] as an indication of a dynam- 
ical phase transition. In this paper we show that 
in driven 3D layered superconductors, the peak 
in differential resistance also coincides with the 
onset of correlation along the c-axis. 

2. MODEL 

Simulations of vortices in 3D layered supercon- 
ductors have been done previously using time- 
dependent Ginzburg-Landau-Lawrence-Doniach 
equations (two layers), fL2) ], the 3D XY model 
flg| l and Langevin dynamics of interacting parti- 
cles JlJJl^]. Here we study the motion of pan- 
cakes vortices in a layered superconductor with 
disorder, with an applied magnetic field in the c- 
direction and with an external homogeneous cur- 
rent in the layers (ab-planes). We use a model in- 
troduced by J.R. Clem for a layered superconduc- 
tor with vortices in the limit of zero Josephson- 
coupling between layers |l6| . We simulate a stack 
of equally spaced superconducting layers with in- 
terlayer periodicity s, each layer containing the 
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same number of pancake vortices. The equa- 
tion of motion for a pancake located in position 
Ri = (ji, Zi) — (%4,yi,ms) (with z-axis in c- 
direction) is: 

where p.y = — | and Zij — \z; t — Zj\ are the 
in-plane and inter-plane distance between pan- 
cakes r ip = |r^ — r p | is the in-plane dis- 
tance between the vortex i and a pinning site 
at R p = (r p ,z.i) (pancakes interact only with 
pinning centers within the same layer), 77 is the 
Bardeen-Stephen friction, and F = -^-J x z is 
the driving force due to an applied in-plane cur- 
rent density J. We consider a random uniform 
distribution of attractive pinning centers in each 
layer with F p = — A p e^^ r ^ r ^ r/r 2 , where r p is 
the pinning range. The magnetic interaction be- 
tween pancakes F v (p, z) = F p (p, z)f is given by: 

F p ( P ,0) = (^/4^ Ap )[i_(A || /A)(l-e-' , / A ")](2) 

F p (p,z) = (^A||/47r 2 A 2 p)[e-*/ A ii - e^] . (3) 

Here, R — \J z 2 + p 2 , X\\ is the penetration length 
parallel to the layers, and A = 2A 2 /s is the 2D 
thin-film screening length. A analogous model to 
Eqs. (2-3) was used in iQ. We normalize length 
scales by Am, energy scales by A v = 0§/47r 2 A, 
and time is normalized by r = i]\ 2 /A v . We con- 
sider N v pancake vortices and N p pinning cen- 
ters per layer in Ni rectangular layers of size 
L x x L y , and the normalized vortex density is 
n v = N v X 2 /L x L y ~ B A 2 /$o- Moving pancake 
vortices induce a total electric field E = — v x z, 
with v = N L N J2i v »- We study the dynamical 
regimes in the velocity-force curve at T = 0, solv- 
ing Eq. (1) for increasing values of F = Fy. We 
consider a constant vortex density n v = 0.1 in 
Ni = 5 layers with L x /L y = a/3/2, s = 0.01, 
and N v = 36 pancake vortices per layer. We 
take a pinning range of r p — 0.2, pinning strengh 
of A p /A v = 0.2, with a density of pinning cen- 
ters n p — 0.65 on each layer. We use periodic 
boundary conditions in all directions and the peri- 
odic long-range in-plane interaction is dealt with 
exactly using an exact and fast converging sum 



|l7|. The equations are integrated with a time 
step of At = O.Olr and averages are evaluated 
in 16384 integration steps after 2000 iterations 
for equilibration (when the total energy reaches a 
stationary value). Each simulation is started at 
F = with an ordered triangular vortex lattice 
(perfectly correlated in c-direction) and slowly in- 
creasing the force in steps of AF = 0.1 up to 
values as high as F = 8. 

3. RESULTS 

We start by looking at the vortex trajectories 
in the steady state phases. In Figure l(a-b) we 
show a top view snapshot of the instantaneous 
pancake configuration for two typical values of 
F. In Figure 2(a-b) we show the vortex trajec- 
tories {Rj(<)} for the same two typical values of 
F by plotting all the positions of the pancakes 
in all the layers for all the time iteration steps. 
In Fig. 3(a) we plot the average vortex velocity 
V = (V y (t)) = (j^Ei^), in the direction of 
the force as a function of F and its corresponding 
derivative dV/dF. We also study the pair distri- 
bution function: 

9(P> n ) = JfJ^C^^p - p i0 )8{ns - Zij )). (4) 

In Fig. 3(b) we plot a correlation parameter along 
c-axis defined as: 

C: = \im g(p,n) , (5) 

as a function of F for n = 1, 2. Below a critical 
force, F c as 0.4, all the pancakes are pinned and 
there is no motion. At the characteristic force, 
F p « 0.8, we observe a peak in the differential re- 
sistance. At F c pancake vortices start to move in 
a few channels , as was also seen in 2D vortex sim- 
ulations 0. A typical situation is shown in Fig. 
2(a). In this plastic flow regime we observe that 
the motion is completely uncorrelated along c- 
direction, with C™ w for F c < F < F p as shown 
in Fig. 3(b). In Fig. 1(a) and 2(a) we see that this 
situation corresponds to a disordered configura- 
tion of pancakes and to an uncorrelated structure 
of plastic channels along c-axis. At F p there is an 
onset of correlation along the c-axis and pancakes 
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vortices start to align forming well defined stacks 
or vortex lines. This onset of c-axis correlation 
corresponds to the transition from plastic flow to 
a moving smectic phase (a complete discussion of 
the translational as well as temporal order will be 
discussed elsewhere 18 1) . For F > F p we observe 
that the structure of smectic channels is very cor- 
related in the c-direction. 
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Figure 2. Pancake trajectories for two typical 
values of F, (a) F c < F = 0.7 < F p , (b) 
< F = 2.5, obtained plotting all the positions 
of the pancakes for all time iteration steps. 



Figure 1. Pancake configuration (top view of the 
layers) for two typical values of F, (a) F c < F = 
0.7 < F p , (b) F p < F = 2.5. Diamonds (O) repre- 
sent pancake vortices and asterisks (*) represent 
pinning centers. 



4. DISCUSSION 

In a system with Ni = 5 layers we have found 
that there is a clear onset of c-axis correlation 
with increasing driving force, where the pancake 
vortices start to align forming well defined stacks 
moving in smectic flow channels. Below this tran- 
sition, in the plastic flow regime, these stacks of 
pancakes are unstable. Also, in jl2j an enhance- 
ment of c-axis correlations with increasing drive 
was observed in a bilayered system. 

We have further found |Q that the in-plane 
properties are in well correspondence with the 
ones obtained in 2D thin films simulations fllQl . A 



better understanding of the effects of c-axis cor- 
relations in pancake motion on each layer can be 
obtained by studying translational and temporal 
order in larger systems, through the analysis of 
the structure factor, voltage noise, and in-plane 
]Tl| ] as well as inter-plane velocity-velocity corre- 
lation functions [Q. 

In conclusion we have analyzed the vortex cor- 
relation along the field direction (c-axis) in the 
velocity-force characteristics at T = and found 
that above the critical current there is an onset 
of c-axis correlation in the transition between a 
plastic flow regime to a smectic flow regime. This 
transition coincides with the peak in the differen- 
tial resistance. Experimentally, this effect could 
be studied with measurements of c-axis resistiv- 
ity as a function of an applied current parallel to 
the layers [fU . 
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Figure 3. (a) Velocity-force curve (voltage- 
current characteristics), left scale, black points, 
dV/dF (differential-resistance), right scale, white 
points, (b) Correlation parameter along c-axis 
C? for n = 1,2. 
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